Analytical results for entanglement in the five-qubit anisotropic Heisenberg model 
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We solve the eigenvalue problem of the five-qubit anisotropic Heisenberg model, without use of 
Bethe's Ansatz, and give analytical results for entanglement and mixedness of two nearest-neighbor 
qubits. The entanglement takes its maximum at A = 1 (A > 1) for the case of zero (finite) 
temperature with A being the anisotropic parameter. In contrast, the mixedness takes its minimum 
at A = 1 (A > 1) for the case of zero (finite) temperature. 
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Recently, the study of entanglement properties of 
many-body systems has received much attention [H|-[2^|. 
To obtain analytical results for entanglement, one may 
consider the case of infinite lattice or a small lattice with 
a few qubits. It is hard to get some analytical results 
between these two extreme cases. 

It was well-known that the anisotropic Heisenberg 
model can be solved formally by Bethe's Ansatz 
method [23l |24| for arbitrary number of qubits N, how- 
ever, we have to solve a set of transcendental equations. 
For N < 7, the isotropic Heisenberg Hamiltonian can be 
analytically solved pa, l2rjj ] . Here, we give the analytical 
results of the eigenvalues of the anisotropic Heisenberg 
model with N = 5, without use of Bethe's Ansatz, from 
which the analytical expressions for entanglement and 
mixedness of two nearest-neighbor qubits are readily ob- 
tained. 

It is interesting to see that the entanglement properties 
of a pair of nearest-neighbor qubits at a finite tempera- 
ture is completely determined by the partition function. 
The entanglement, quantified by the concurrence [27| . 
relates to the partition function Z via [28l [29l l30| 



C = max 



U 
N 



AG, 



with 



U = — 



dlnZ 



dp 



(1) 



(2) 



being the internal energy, and 

G zz = Tr[eM~PH)a lz a l+lz ]/Z = (3) 

being the correlation function. Here, (3 — 1/T and the 
Boltzmann's constant k = 1. Thus, once we know the 
eigenenergies versus the temperature and the anisotropic 
parameter, we can completely determine the entangle- 
ment. 

There exists another concept, the mixedness of a state, 
is central in quantum information theory |3l| . For in- 
stance, Bose and Vedral have shown that entangled states 
become useless for quantum teleportation on exceeding 
a certain degree of mixedness [33 . Mixedness is also re- 
lated to quantum entanglement. We will study both the 
entanglement and mixedness properties. 



Eigenvalue problem. The anisotropic Heisenberg 
Hamiltonian is given by 



J 



N 



H =— ^2 (1 + Vix^i+ix + o'lyO'i+iy + Aa lz a l+ i z ) 



A- 1 



(4) 



where Sjj+i —\ (1 + <7j • 0j+i) is the swap operator be- 
tween qubit i and j , <Xj = (cti X , o~i y , <7j z ) is the vector of 
Pauli matrices, and J is the exchange constant. We have 
assumed the periodic boundary condition, i.e., N+l = 1. 
In the following discussions, we also assume J = 1 (anti- 
ferromagnetic case) and A > 0. 

Since we impose the periodic boundary condition, the 
Hamiltonian is translational invariant, i.e., [H, T] = 0, 
where T is the cyclic right shift operator defined as 



T\mt, ■ ■ ■ ,mtf-i,m N ) = \m Nl mi 



,mjv_i). (5) 



The translational invariant symmetry can be used to re- 
duce the Hamiltonian matrix to smaller submatrices by 
a factor of H^. 

Now we focus our attention to five-qubit settings, and 
solve the eigenvalue problem of the anisotropic Heisen- 
berg model. Since [H, J z ] = 0, the whole 32-dimentional 
Hilbert space can be divided into invariant subspaces 
spanned by vectors with a fixed number of reversed spins. 
Then, the largest subspace is 10-dimensional with 2 or 3 
reversed spins. Here, J z = Y^=i °~iz/2- Dire to the sym- 
metry [H, iy = 0, it is sufficient to solve the eigenvalue 
problems in the subspaces with r reversed spins, where 
r e {0, 1, 2} and T, x — a x ® a x ® a x ® a x <£> o x . By us- 
ing the translational invariance, we can further reduce 
the Hamiltonian matrix to 2 x 2 submatrices, and the 
eigenvalue problem can be readily solved. 

The subspace with r = only contains one vector 
1 00000), which is the eigenvector with eigenvalue 



E = 5[l + (A-l)/2]. 



(6) 



The subspace with r = 1 is spanned by five basis vectors 
{T"|10000),n G {0,1,2,3,4}}. Considering the transla- 
tional invariance of the Hamiltonian, we choose another 



2 



basis given by 



4 

E 

n=0 



<T™ 1 10000), 



(7) 



where u>k = e 5- ^,^ S {1,2,3,4,5}. It can be checked 
that states \(pk) are eigenstates of T with eigenvalues oj^ 1 , 
and are also eigenstates of Hamiltonian H with eigenval- 
ues given by 



Ei,. 



A- 1 



2cos(2fc7r/5). 



(8) 



For the case of r = 2, we choose the following basis for 
10-dimensional subspace 



\>l>h 



J2^T n \ 11000), 

n=0 
4 

£ wjJT" 1 10100). 



(9) 



n=0 



States j^fe) and \<j>k) span an invariant 2x2 subspace 
under the action of Hamiltonian H . In this subspace, 
the Hamiltonian can be written as 



Hi. — 



2 ± 
U k 1 



3(A-1) 



(10) 



Then, from the above equation, the eigenvalues are ob- 
tained as 



E- 



2cos(4/br/5) 



2.fe± = 

z 

±a/[A - cos(4/br/5)] 2 + 2[1 + cos(2/br/5)]. (11) 

Thus, all eigenvalues are obtained for the five-spin 
anisotropic Heisenberg models. We see that eigenstates 
are at least two-fold degenerate due to the symmetry 
[if, H x ] = 0. Although the eigenstates can be easily ob- 
tained, they are not given explicitly here as the knowl- 
edge of eigenvalues is sufficient for discussions of entan- 
glement and mixedness properties. 

Entanglement and mixedness. From Eqs. 0-©, the 
concurrence is determined by the partition function and 
its derivative with respect to the anisotropic parame- 
ter A. As we have obtained all eigenvalues of five-spin 
Hamiltonian if, from Eqs. ©-©i it follows that 



Z 



-0E O 



E (*- 0El ' k 

fe=l 



-0E 2 



-0E 2 



(12) 

Substituting the above equation into Eqs. 10) and Q 
yields 



Eoe-^ 
E2,k+e 



5 

< E 

fc=i 



(-Bi, fee" 



E 2 ,k-e 



-0Ei,k 



(13) 
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FIG. 1: The concurrence versus A for different temperatures. 



and 



Gzz ~hZ 



-0E O 



E (r pEl ' k 

k=l 



dE 2: k+ 
<9A 



I3E 2 



with 



dE 



2,fc± 



dA 



1 



dA 



A - cos(4/c7r/5) 



(14) 



cos(4/c7r/5)] 2 + 2[1 



cos(2A:7r/5)] 

Then, the analytical expressions of the internal engery U 
and the correlation function G zz are obtained, and thus, 
according to Eq. QJ, we get the analytical expression of 
the concurrence for the thermal state. 

At zero temperature, the system is in the ground state, 
and for this case, Eq. (JJJ reduces to 



C 



0. 



,E, 



N 



A8E gs 
N dA 



1 dE gs 
N dA 



(16) 

where E gs denotes the ground-state energy and N = 5. 
In the derivation of the above equation, we have used the 
relation 



G zz = — 



N K dA' 



2 dE a 



N dA 



(17) 



For our five-qubit 



which is valid for the ground state, 
model, E gs = E 2 ,i- 

spect to A is given by Eq. (TSJ. Thus, the analytical 
expression for the entanglement of ground state is ob- 



and the derivative of E gs with re- 



tained, 
to 



If 



N 



A 9E g 
N dA 



< 0, then Eq. I|16|) reduces 



C 



E. 



N 



ldE a 



N dA 



(18) 



where we have ignored the max function. Taking deriva- 
tive with respect to A on both sides of the above equation 
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FIG. 2: The linear entropy versus A for different tempera- 
tures. 



FIG. 3: The concurrence (square and diamond lines) and lin- 
ear entropy (circle and star lines) versus A for different tem- 
peratures in the four-qubit model. 



leads to 



dC_ 
OA 



A - 1 d 2 E gs 
N d 2 A 



(19) 



Then, it is direct to check from the above equation and 
the ground state energy E 2 .\- that the derivative 
is zero when A = 1. Thus, the concurrence takes its 
extreme value at the point of A = 1 for the ground state. 

From the analytical results for the concurrence, we nu- 
merically plot the concurrence versus the anisotropic pa- 
rameter for different temperatures in Fig. 1. We observe 
that the concurrence takes its maximum when A = 1. 
This point correspond to the critical point of metal- 
insulation transition p9l |. However, for finite tempera- 
ture, the concurrence reaches its maximum when A > 1. 
For finite temperatures (for instance, T — 2.0), we find a 
threshold value of the anisotropic parameter A t h, before 
which there is no pairwise entanglement. The threshold 
value increase as temperature increases. 

Next, we study mixedness properties of the thermal 
state and ground state. The mixedness of a state g can be 
quantified by the linear entropy given by El = 1 — Tr(g 2 ). 
Then, for arbitrary number of qubits, the linear entropy 
of the state of two nearest qubits is given by 



E L {T) = l-~ 



2(--- 
\N 2 



AG, 



for the thermal state, and 



E L (T = 0)=1-- 



Egs_ 

N 

4 /dE g3 

N 2 V dA 



G z 



A BE, 



1 



(20) 



N dA 



1 



(21) 



for the ground state, respectively. We see that the mixed- 
ness of the thermal state is also completely determined by 
the partition function, and the mixedness of the ground 
state is determined by the ground-state energy and its 



first-order derivative with respect to the anisotropic pa- 
rameter. Thus, the analytical expressions of the linear 
entropy are obtained. 

In Fig. 2, we numerically calculated the linear entropy 
versus A for different temperatures. In contrast to the 
entanglement, the mixedness of the ground state displays 
a minimum when A = 1. For the case of finite temper- 
atures, the mixedness takes its minimum when A > 1. 
For lower temperatures (for instance, T < 0.2), numer- 
ical results show that the maximum of the concurrence 
occurs nearly at the same value of A as the minimum of 
the mixedness. It seems that the more the pairwise en- 
tanglement, the less the mixedness. However, for higher 
temperatures, the maximum of the concurrence and the 
minimum of the mixedness do not occur at the same A. 
For instance, when T = 1.5, the concurrence takes its 
maximum at A = 3.1037, and the mixedness takes its 
minimum at A = 3.8525. This signifies that it is not al- 
ways true that the more the pairwise entanglement and 
the less the mixedness. 

For the case of four qubits, the anisotropic Heisenberg 
model is also exactly solvable using the same method as 
above. Here, we make a comparison of the four-qubit and 
five-qubit cases. The exact ground-state energy and its 
derivative with respect to A are given by 



2- A- VA 2 



T 9E gs 
8 ' ~JK 



1 



VA 2 



(22) 



Substituting the above equation into Eqs. I|16fl and 121(1 . 
we obtain the ground-state concurrence and linear en- 
tropy as 



C(T 
E L (T 



0) 



A- 



4VA 5 



A 



_11 

) ~16 ~ 8^ 



32 



~8 16(A 2 



(23) 
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From the above analytical expressions, it is straightfor- 
ward to check that the concurrence takes its maximum 
and the linear entropy takes its minimum at the point 
of A = 1, which is exactly the same feature as that in 
the five-qubit model. For finite temperatures, in Fig. 3, 
we give numerical calculations of the concurrence and 
the linear entropy. We see that they displays similar be- 
haviours as those in the five-qubit model. For instance, 
for T > 0, the maximum pairwise entanglement and the 
minimum mixedness occur at A > 1. 

Conclusion. In conclusion, we have obtained the an- 
alytical results for the entanglement and mixedness in 
the five-qubit anisotropic Heisenberg model. The ex- 
act eigenspectrum is obtained, and entanglement and 
mixedness properties can be completely determined by 
the eigenvalues of the system, irrespective of the eigen- 
states. The method adopted here can be applied to the 
anisotropic Heisenberg model with more than five qubits 
(for instance, 6 or 7 qubits). However, the analytical ex- 
pressions for eigenvalues, entanglement, and mixedness 
are expected to be more complicated. 

We have made numerical calculations, and show that 
the entanglement takes its maximum at A = 1 (A > 1) 
for the case of zero (finite) temperature. In contrast, 



the mixedness takes its minimum at A = 1 (A > 1) for 
the case of zero (finite) temperature. From our analy- 
sis, we conjecture that it is a general feature that at zero 
temperature the entanglement takes its maximum and 
the mixedness takes its minimum when A = 1 for any 
number of qubits. This conjecture is supported by our 
analytical results for four and five qubits, and by numer- 
ical results for the number of qubits being as large as 
1280 |2^|. The Heisenberg chains not only displays rich 
entanglement features, but also have useful applications 
such as the quantum state transfer |34[ . Experimentally, 
it was found that entanglement is crucial to describing 
magnetic behaviors in a quantum spin system |35|. So, 
the study of entanglement and mixedness properties in 
the Heisenberg models will strength our understanding 
of other quantum features of magnetic systems. 
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